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ABSTRACT 

A five-parameter fitting formula for the line-of-sight stellar velocity distributions of 
steady state systems is proposed. It can faithfully reproduce velocity distributions of 
theoretical models ranging from nearly Gaussian profiles to strongly skewed or mildly 
double-peaked profiles. In contrast to van der Marel and Eranx (1993) and Kuijken and 
Merrifield (1993), the line profiles are required to have neither multi-peaks nor negative 
velocity wings and most informations can be recovered from no more than five physically 
meaningful fitting parameters. 
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1. Introduction 

Exactly Gaussian stellar velocity distributions exist only in simplified theoretical 
models. Line-of-siglit velocity distributions (LOSVDs) of realistic stellar systems, which 
are implied by broadened spectral features in a galaxy spectrum by the Doppler effect, 
commonly show deviations from Gaussian in terms of a non-zero skewness and kurtosis, 
and a fair fraction of them can be slightly double-peaked. When fitting broadened spectral 
lines from observation, it is important to parametrize the LOSVD more generally than 
a three-parameter Gaussian because the deviations from Gaussian carry as much critical 
information of the stellar dynamics of the underlying system as the line-strength, the mean 
velocity and the dispersion. For example, LOSVDs have been used to identify kinematically 
distinct cores in some elliptical galaxies (Franx and lUingworth 1988, Rix and White 1992) 
and set non-trivial limit on the dark halo in an FO galaxy (CaroUo et al. 1995). Similarly in 
two spiral galaxies analysis of LOSVD have revealed that a fraction of stars in the disc are 
counter-rotating (Rix et al. 1992, Merrifield and Kuijken 1994). Another recent application 
is the detection of a bar in two edge-on spiral galaxies with boxy bulges based on the 
presence of double-peaked stellar LOSVD (Kuijken and Merrifield 1995). 

When fitting a typical galaxy spectrum in the presence of noise and systematic fitting 
errors, it is likely that the well-determined parameters is only a handful. It appears to be 
more sensible to fit the observation with a simplely parametrized velocity profile than a 
profile with many parameters or a fully non-parametric profile. If all realistic profiles were 
exactly Gaussian, then the number of meaningful parameters would be only three, namely, 
the amplitude, the mean velocity and the dispersion. If realistic profiles were genericly very 
close to Gaussian, then the number of meaningful parameters would be no more than five; 
the two additional parameters would be the skewness and the kurtosis. In this paper we 
show that even for profiles with large deviation from Gaussian can still be specified by no 
more than five parameters. 

Realistic profiles are likely much more restricted than an arbitrary positive function 
of velocity, because the velocity and density distributions of a steady state system satisfy 
several equations including the Vlasov equation and the Poisson equation. Previous methods 
are often designed to cover all smooth positive and negative mathematical functions in a 
complete basis set as the Gauss-Hermite expansion method by Gerhard (1993) and van der 
Marel and Franx (1993, hereafter vdMF), or smooth positive functions as the unresolved 
Gaussian decomposition method by Kuijken and Merrifield (1993, hereafter KM). These 
parametrizations also tolerate multiple peaks, peaks at the wings and even negative wings 
in the profile, which are likely spurious due to noise and artifacts of the parametrizations. 
Many free parameters are often invoked to fit theoretical profiles or observed spectrum. 
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typically about 5 to 9 parameters for vdMF, and on the order of 20 parameters for KM. 

Different from these methods, our philosophy is to maximally exploit the a priori 
information of realistic profiles so as to represent a general but realistic profile with only 
a few parameters. The parametrization is tailored to model the majority of realistic line 
profiles which have the following properties by assumption. 

(I) They are positive and smooth everywhere in velocity. 

(II) They have no more than three points with vanishing first derivative excluding 
at infinity, namely, they are either single-peaked profiles with up to two fiat shoulders or 
double-peaked profiles. 

(III) The deviation from a Gaussian at the wings is much smaller than the peak 
amplitude and/or is not significantly above the noise level. 

These conditions are fairly general, and are often met by line profiles of theoretical 
models and of observed systems (see references in the beginning of the Introduction) where 
any strong deviation from Gaussian is near the systematic velocity. They also include most 
of the interesting classes of line profiles from slightly skewed profiles to mildly double-peaked 
profiles. 

In this paper we give parametrizations with the above properties, and compare with 
previous methods. Both Gerhard (1993) and vdMF use the Gaussian-Hermite expansion 
parametrization; the former uses a fixed Gaussian as the zero order term while the latter 
uses the best fitting Gaussian. Our comparison will be mostly with the more related vdMF 
method. 

The paper is organized as following. §2 gives the parametrization and its properties. It 
is then tested with synthesized profiles from theoretical models in §3. The conclusion and a 
brief discussion of possible generalizations are given in §4. 



Motivated by the fact that Gaussian is still a fair approximation to realistic profiles, 
we write a general profile L{v) as the following. 



2. 



The Parametrization 



2.1. A modified Gaussian-Hermite parametrization 
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where 

v-V 

,„ = — , (2 

is the rescaled velocity, 7, a are the hne strength, mean velocity and dispersion of the 
best fitting Gaussian. Q{w) specifies the amount of deviation from the Gaussian. 

In vdMF Q{w) is given by an expansion on the complete Gaussian-Hermite function 
set Hn{w) with n > 3 (see their equation 9). Namely, 

N 



QvdMpiw) = J2^nHn{w), (3) 

n=3 

where c„ is a set of fitting coefficients for the amplitude of deviations from Gaussian, and 
Hn{w) is the n-th Gaussian-Hermite function defined in vdMF, which satisfies 

Hm{y)Hn{y)e ^ dy = VT^^mn , (4) 



and are more specificly given by 

Hsiy) = ^iy'-ly), 

and 

H,iy) = -^iy'-2y' + j). (6) 

QvdMF{w) is essentially a polynomial function of w of order three or higher, which like 
virtually all polynomial functions is wildly oscillating and asymptoticly reaches positive 
or negative infinity. An example is shown by the dashed lines (labeled 6 = 0) in Figure I 
for (up to a constant factor) H4[w) (right panel) and H^{w) (left panel). Due to these 
behaviours of Q{w), the line profile L{v) often has multiple peaks or negative wings or 
artificial peaks at the wings. 

For realistic profiles with the three properties given in the Introduction, we propose 
to modify the vdMF parametrization and to regularize Q{w) and L{v) at the wings by 
introducing Gaussian damping terms exp( — 6„t/;^/2) with the constants 6„ > 0. Generally 
we write 

Q{w) = ^Cnane~ 2 Hn{anw), a„ = yl + Y. (7) 

This parametrization is designed so that it reduces to QvdMF{w) if 6„ = for n = l^N. 
In practice we find it is sufficient to truncate the expansion at = 4. As a result the 
total number of fitting parameters are only five, namely, (7, cr, C3, C4). We further set the 
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constants 63 = 64 = 6 and impose one of the three cases of damping at the wings: zero 
damping [b = 0), medium damping [b = 1) and strong damping [b = 2); the zero damping 
corresponds to vdMF and merely serves as a reference case. 

Specificly the proposed parametrization of hne profiles is as follows, 

L{v) = -^e-^[l + Q{w)] (8) 

e 2 1^ + ae 2 ['j:iH3{aw) + -f4H4{aw)\} (9) 



where 



b V - V 

fee {0,1,2}, a = Wl + -, 73 = C37, 74 = C47, w = ——. (10) 



2.2. Properties similar to Gaussian-Hermite expansions 



Similar to vdMF, the parameters (7, cr, C3, C4) used here are distinct moments of the 
line profile L[v). They are related to L[v) by the following equations. 

r- /■ + °° »2 

V2 / L(v)e-—dv, (11) 



7 

— 00 

1 /■ + °° u,2 

V = - vL{v)e-—dv, (12) 

7 

1 ft -|- Q)Q 2 

a = [- (v -VYL{v)e-'^dv]2, (13) 
7 J-00 

\/2 / + °° (l + b)ti,2 

C3 = — — / L{v)e 2 H'i{aw)dv^ (14) 
7K3 J-00 

-^2 /■ + °° (l + b)ti.2 

C4 = — — / L{v)e 2 HAaw)dv . (15) 
7K4 J-00 

The constants and A;4 are given in the Appendix. These relations can be easily shown 
if one applies equation (4) and also notes that cross terms between an even and an odd 
Gaussian-Hermite function are always zero and that by design 

H^{w)Q{w)e-'" dw = {)ioi m = {),l,2. (16) 



The above equations imply that our parametrization as in equations (9)-(10) preserves 
the nice properties of vdMF. First (7, cr) are the parameters for the best fitting Gaussian. 
The parameters C3 and C4 describe the anti-symmetric and the symmetric deviations from 
Gaussian; they reduce to and /i4 if 6 = or if the deviation from Gaussian is small. 
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Second the basis set in our parametrization is an orthogonal basis set and the 
parameters correlate only weakly (due to non-linearity). Fixing V and a so that the 
parametrization is a linear function of the fitting parameters (7,73,74), the correlation 
matrix of (7,73,74) is diagonal and equals to the identity matrix. The nearly diagonal 
covariance matrix simplifies the error analysis. 



Third it is also easy to implement the convolution with template spectrum with our 
parametrization. The Fourier transformation of L[v) is simplely 

''^^ Liv^e'^^dv = 76"^^' + e"^^'{-z739"2[i73(-/) + VSr^] (17) 

-00 (I 



a a V Z 



where 



b b 

q=l+r = 1 + -, a=\ l + -. 18 

^ 2 + 6' V 2 



2.3. New properties 



The main benefit of the modification is that the spurious positive or negative peaks at 
the wings in vdMF can be effectively suppressed in this parametrization. With damping 
Q{w) is finite everywhere and asymptoticly approaches zero because mathematically the 
Gaussian damping term dominate the polynomials at large velocities. Hence for sufficiently 
small c„, l-\-Q{w) is positive everywhere and unity at infinity. The resulting profile L{v) can 
deviate from Gaussian only near the systematic velocity, not at the wings. The behaviours 
of Q{w) in two cases 6=1 and 6 = 2 as well as the 6 = reference case are illustrated by the 
solid lines in Figure 1, which show (up to a constant factor) Hn{w^Jl + b/2) exp[ — bw^ /2) 
versus w for n = 3 (right panel) and n = 4 (left panel). 

Fig. 2 shows some typical profiles in each parametrization. While vdMF (as shown by 
the 6 = undamped case) genericly contains negative wings which is significant even for 
slightly fiat-topped profiles or skewed profiles, the damped cases (6 = 1 or 6 = 2) show no 
negative wings. The 6=1 profiles range from nearly Gaussian to mildly double-peaked 
profiles, including fiat-topped and strongly skewed profiles. The 6 = 2 profiles can also be 
strongly double-peaked. Still we prefer the parametrization for the 6=1 case with equation 
(9)-(10) because strongly two-peaked profiles are relatively rare compared to skewed and 
mildly double-peaked profiles. In general the profiles with 6=1 also approaches to Gaussian 
at the wings slower than the 6 = 2 profiles, hence can incorporate non-Gaussian wings, e.g., 
a power-law fall-off, near 2.5(7. 
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The shape of our profiles is completely specified by C3 and C4. Fig. 3 and Fig. 4 show the 
regions of C3 — C4 that prescribe smooth positive profiles with one, or two, or three peaks for 
the 6=1 case and for the b = 2 case. When C3 and C4 is small (approximately < C4 < 0.3 
and |c3 lb C4I < 0.3), the profile is single peaked and (7, cr, C3, C4) are approximately the 
amplitude, the mean velocity, the dispersion, and up to a constant (see Appendix) the 
skewness and the kurtosis of the profile. For sufficiently big C3 or C4, the profile will contain 
negative wings or multiple (> 3) peaks. For intermediate C3 and C4, the profile becomes 
double-peaked. The meaning of the parameters are less direct than, e.g., the ratio and the 
offset between the two peaks. In these cases one can use Table I (if 6 = 1) and Table II (if 
6 = 2) to convert C3 and C4 to more direct characteristic shape parameters of the profile and 
scale the physical dimensions with (7, cr). 

Comparing with previous parametrizations, our parametrization as in equations 
(9)-(10) not only can model a wide range of realistic profiles, but also has many nice 
properties, including using the least number of free parameters. All information of the 
line profile is contained in five numbers. Because damping makes it possible to fit a range 
of realistic profiles without introducing spurious wings, there is little need to include 
higher order Gaussian-Hermite terms such as H^{w) or He[w) in our parametrization. 
By contrast, one often needs to keep higher order expansion terms Hn[w) in vdMF 
method with + 1 = 5 to 9 parameters to hide the negative wings effectively and to 
fit double-peaked profiles. The multi-Gaussian decompositions (KM) and non-parametric 
fitting effectively invoke even more parameters, on the order of the number of velocity bins. 
The double-Gaussian decomposition covers similar range of profiles as our parametrization, 
but the former invokes six parameters, which often do not have clear physical meanings and 
are often strongly correlated. 



3. Test with synthesized profiles 



In this section we apply the parametrization to some synthesized profiles from simple 
theoretical models of spheroid and disc systems. 

Figure 5 and Figure 6 show a few line profiles of some analytical spheroidal models 
inside a fiattened logarithmic dark halo potential by Fvans (1994). These models are 
fiattened and have power-law density profiles and anisotropic velocity ellipsoids. When 
viewed edge-on on the major axis, the projected line-of-sight velocity distribution of a 
non-rotating anisotropic model can be written analytically as following. 

Liv) = e-^+Se-^^\l + kv'), S = '-^{^^^^K ^ = (p " l)(p + 2)-^, (19) 
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where the dark halo core radius and terminal velocity have been set to unity, p is the 
power-law slope of the stellar spheroid, s in Evans's notation or 6 are the anisotropy 
constant, R is the major axis distance. Since a rigorous expression for rotating systems is 
complicated, we simulate effects of rotation qualitatively by replacing the first term 
with e 2 ^ where t is a rotation measure. 

We have computed a few profiles at different locations [R = f,2,oo) in the analytical 
models with different amount of anisotropy [6 = 0.2, f, 2), rotation [t = 0,0.2,0.4) and 
power-law slope [p = 3 for Figure 5 and p = A for Figure 6) based on equation (f9). These 
profiles are then fitted using both the vdMF method (the heavy dashed line) and our 
method with b = 1 (the solid line). We find that our parametrization fits the analytical 
profiles better than vdMF in all above cases although larger residuals of our fit start to 
show up outside the core of the steeper power-law case (Figure 6). 

The analytical profiles cover a wide range in shape. First the deviation from Gaussian 
becomes smaller closer to the core. Second at fixed position R the profiles change 
continuously from a double-peak structure to a nearly Gaussian shape for decreasing values 
of anisotropy 6 and/or power-law slope p. These profiles are similar in variety to those of 
other analytical spheroidal models, e.g., those shown in Gerhard (f993). It is also likely 
that these profiles span the range of profiles of real bulges or ellipticals, which are hot and 
mildly fiattened and anisotropic systems with typical power-law slope less than or equal to 
4. While some line profiles of systems with very big anisotropy, fiattening and power-law 
slope are not well-fitted by our parametrization as shown in Figure 6, these systems may be 
intrinsicly rare in nature due to formation processes and constraints from equilibrium and 
stability. 

It is necessary to caution that the two-integral models used as test cases here are 
likely too simple for real elliptical galaxies, where the third integral always plays a role. 
Based on the observed line profiles of a sample of 50 or so ellipticals Bender, Saglia and 
Gerhard (1994) found that they deviate from a Gaussian by no more than 10%, which 
would argue that real profiles are typically single-peaked. Deviations are more prominent in 
the asymmetric part than the symmetric part with |c3| j/is] < 0.15 and |c4| \h4\ < 0.05. 
Profiles with a strong fiat-top (c4 /i4 < —0.05) are absent in their sample, which they take 
as evidences against two-integral fiattened non-rotating models being the right dynamical 
model. 

We have also tested the parametrization with simple models of disc systems. We 
synthesize line profiles of a disc plus bulge systems by superimposing a wider Gaussian 
with a shifted narrow Gaussian. Similarly we synthesize line profiles of counter-rotating 
disc systems by superimposing two well-separated Gaussian profiles with equal dispersion. 
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In both cases we vary the ratio of the two Gaussians to simulate profiles at different major 
axis distance. Figure 7 shows the fit to these double-Gaussian line profiles. The results are 
consistent with the finding with analytical spheroid models. Except for some extreme cases 
as shown in panel 4 where the amplitude ratio is very small, most line profiles are well fit 
with our parametrization. 



4. Conclusion and Possible Generalizations 



In summary, we describe a method to derive line profiles without spurious wings. 
The method is based on modifying the Gaussian-Hermite expansion method by Gerhard 
(1993) and vdMF. The main application of our method is to recover the commonly seen 
strongly skewed, fiat-topped or weakly double-peaked profiles. The power of the proposed 
parametrization remains to be tested in fitting observed spectrum. 

It is possible to generalize the proposed parametrization. In particular the lowest order 
term needs not be a Gaussian. Many symmetric functions which fall off steeply at large 
velocities would also serve. The following parametrization, for example, is found to work 
very well. 

L{v) = sech—[l + {siw' + S2w'^)] (20) 
where Si and S2 are fitting coefficients, and 

w = , w = ^. zi 

'2 

Fig. 8 shows fits to a theoretical line profile by five different parametrizations, and their 
residuals. The theoretical profile is the symmetric double peaked profile shown in Figure 3d 
of Fvans (1994), which is the profile at one core radius for a non-rotating Fvans model with 
the power-law slope p = = 3.5 and 6 = 0.9. Clearly the vdMF method including H4[y) 
(the dashed line) or both H4[y) and He[y) (the heavy dashed line) cannot fit the double 
peaks and have significant negative wings. The rest three parametrizations fit roughly 
equally well with the double-Gaussian (the dotted line) and the Sech-parametrization (the 
faint solid line) being even better than our preferred parametrization (the heavy solid 
line) in this case. Also curiously our parametrization recovers the true dispersion and 
kurtosis of the symmetric line profile better than both the double-Gaussian parametrization 
and the Sech-parametrization as shown in Table 3; the truncated vdMF expansions are 
also disfavored in this criteria. In general we perfer the parametrization in equation 
(9)-(I0) because it is easy to implement and because the involved parameters are few and 
uncorrelated. 
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5. Appendix 

In this appendix we give detailed expressions lor a lew quantities used in the main text. 
(1) The constants and k4 in equations (14) to (15) are given by the following, 

1 r+oo 

ks = ^ e-^^^ [H3{aw)]\dw (1) 

'it J-oo 



(1 + b)-'^^{l + b/2f^{l + b + 56V8) (2) 
1 lor 6 = 0, (3) 

-Valor 6 = 1, (4) 



256 

iiy6 1or6 = 2, (5) 
k^ = ^ e-^^^ [H4aw)]\dw (6) 

'it J-oo 



(1 + 6)-^/2(l + b/2fil + 2b + 9b' /i + 56^4 + 356^128) (7) 

1 lor 6 = 0, (8) 
8fi7 

lor 6=1, (9) 

4096 ' ^ ' 

227 

y6 1or6 = 2, (10) 



1944 



where a = y 1 + 6/2. 

(2) The following are the asymptotic expressions lor the line strength, the mean 
velocity, the dispersion, the skewness and the kurtosis ol our parametrized profile in the 
limit that the deviation Irom Gaussian is small. 

r+oo 

Line strength 7 = / L{v)dv (11) 



(1 + 6/2)1/2 

= + Vt^^^^] (12) 

= 7(1 + iy6c4) lor6 = (13) 

= 7(1 + A^C4) lor6= 1 (14) 

= 7(l + ^C4)lor6 = 2 (15) 
io 

r+oo 

Mean velocity y = / vL{v)/^dv (16) 
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Table 1. Characteristic parameters of a single or double peaked profile on the grid of 

C3 — C4 in the 6=1 case 
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Note. — For each pair of C3,C4, the table gives the peak ratios and offset 
between the second peak with the main peak and ), the dip with the 

mam peak {§- and -^^^^), the main peak with the best fitting Gaussian (Pi, 
Vi), the mean velocity, the dispersion, the skewness and the kurtosis of the 
profile. 
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Table 2. Same as Table I, but for 6 = 2 
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0.54 , 
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0.54 , 


0.08 


1 


16 
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04 


1 


01 


0.03 , -0.05 





20 


0.30 


0.46 , 


1.28 


0.38 , 


0.44 


1 


30 


1 


16 





04 


1 


04 


0.02 , -0.11 





40 


-0.10 


0.59 , 


1.32 


0.46 , 


0.52 


1 
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98 
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1 


30 


1 
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0.56 


1 


40 


1 
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09 


1 


03 
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60 


0.10 


0.52 , 


1.44 


0.25 , 


0.64 


1 


46 


1 


12 





13 


1 


00 


0.08 , -0.09 





60 


0.30 


0.50 , 


1.44 


0.18 , 


0.64 


1 


53 


1 


20 





13 


1 


03 


0.06 , -0.14 



Table 3. Fitted parameters for a symmetric line profile of Evans model by different 

methods. 



Dispersion Kurtosis 



Evans model 


0.60 


-0.14 


h = 1 


0.62 


-0.10 


Sech 


0.59 


-0.50 


2-Gaussian 


0.58 


-0.59 


vdME, N = 4 


0.45 


-11.78 


vdME, N = 6 


0.65 


1.94 



Note. — compares the dispersion and kurtosis of a line profile from the 
Evans model with those derived from different methods. 
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Fig. 1. — The two panels show the fourth order symmetric and the third order anti-symmetric 
Gauss-Hermite polynomials with different amount of damping b = 0,1,2. The dashed line 
(labeled 6 = 0) corresponds to the no-damping case as in vdMF. More specificly the curves 
show (up to a constant factor) Hn{w^l -\- b/2) ex])[ — bw^ /2) versus w for n = 4 (left panel) 
and n = 3 (right panel). 
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Fig. 2. — The upper, the middle and the lower panels show some typical symmetric and 
asymmetric profiles for b = 0,1,2. The solid lines are double-peaked profiles, the dashed 
lines are strongly squared or sharply peaked or skewed profiles. 



- 16 - 





o single peaked 






+ double peaked 






A triple peaked 






AAA 






AAAAAAA 






AAAAAAAAAAA 






^AAAAA_AAAAAA 






AAAAAAAAAAAAAAA 






AAAAAAAAAAA AAAAAA 






AA++++AAA++++AA 






AA+++++0+++++AA 






A++++000++++A 






A++00000++A 






++00000++ 






AOO+OOA 






+ + + + + 
+ 












-2-10 1 2 

c3 



Fig. 3. — shows the regions of C3 — C4 plane which would yield positive profiles of one, two 
and three peaks. The damping parameter 6 = 1. 
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Fig. 4. — same as the previous figure, but for tlie damping parameter 6 = 2. 



Fig. 5. — shows the fits to fine profiles of a simple power-law galaxy model by Evans (1994) 
with the power-law slope p = 3 (the plus symbols) using our method with 6=1 (solid lines) 
and vdMF method with 6 = (dashed lines). The profiles are at different radius R and with 
different amount of rotation t and anisotropy 6. 




Fig. 6. — same as the previous figure, but for tlie power-iaw siope p = 4. 
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Fig. 7. — shows the fits to double- Gaussian fine profiles (the plus symbols) using our method 
with 6=1 (solid lines) and vdMF method with 6 = (dashed lines). Each panel shows one 
case with two counter-rotating Gaussian discs (cruise = 1 and Vrot = il.35) and one case 
with the direct disc plus a non-rotating Gaussian bulge with (Tbuige = 2. The percentage of 
either the counter-rotating disc or the non-rotating bulge is as labeled. 




Fig. 8. — shows the fits to a symmetric double-peaked profile from Evans' model (the 
plus symbols) by five different parametrizations. The heavy solid line and the solid line 
show our preferred parametrization as in equations (9)-(f0) with b = 1 and an alternative 
parametrization as in equation (20). The dashed line and the heavy dashed line are fits with 
vdMF method as in equations (f), (2) and (3) including up to H4[y) and He[y) respectively. 
The dotted line shows a fit by double- Gaussian parametrization. 



